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Abstract 

Unresolved coalescing binary neutron stars in the Galaxy and beyond produce a 
stochastic gravitational wave background within the LISA frequency band, which can 
be potentially dangerous for possible detection of fundamentally more interesting relic 
On ! cosmological backgrounds. Here we address the question what effects the unavoidable 

eccentricity of orbits of these systems should have on the properties of this backgrounds. 
In particularly, we examine starting from which frequency one-year observations with 
ON , LISA can be secure from contamination by the noise produced by galactic binary 

| neutron stars in eccentric orbits. We come to the conclusion that harmonics from such 

■ binaries do not contribute significantly above ~ 10 -3 Hz. 

o : 

1 Introduction 

Cu- 
lt is generally recognized that in the 21st century gravitational wave astronomy will open up a 

new window to study physical processes in strong gravitational field and in the early Universe 

(see Grishchuk et al. 2001 for a recent review). Among the potential sources of gravitational 

waves (GW) for the forthcoming GW detectors, relic stochastic GW backgrounds, coming up 

. £h [ from the early Universe, play a special role due to their fundamental importance and potential 

detectability by the space-born laser interferometer LISA (Grishchuk 1997, Grishchuk et al. 

S ; 2001). 

The existence of relic gravitational waves is a consequence of quite general assumptions. 
In the framework of the conventional cosmological picture, the strong variable gravitational 
field of the early Universe amplifies the inevitable zero-point quantum oscillations of the 
gravitational waves and produce a stochastic background of relic gravitational waves in a 
wide frequency band from 10~ 18 Hz to 10 10 Hz measurable today. The detection of relic 
gravitational waves is a feasible way to learn about the evolution of the very early Universe, 
up to the limit of Planck era and Big Bang. 

The expected level of this background is model-dependent so its robust estimation is 
complicated. Nevertheless, if the level of some background turns out to be higher than the 
detector sensitivity limit in a given frequency band, it can be directly detected using only one 
interferometer LISA, provided that other backgrounds do not contribute at these frequencies 
(for more deep discussion of the LISA capabilities see recent papers by Cornish and Larson 
(2001), Cornish (2001a,b), and Ungarelli and Vecchio (2001)). 
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Among various astrophysical GW backgrounds, those produced by several classes of merg- 
ing binary compact stars (white dwarfs, neutron stars, black holes and their combinations) 
dominate in the LISA frequency band 10 -5 — 10 -2 Hz (e.g. Hils and Bender 1998, Kosenko 
and Postnov 1999, Ferrari et al. 2001). A feature of all such backgrounds is that their level 
(in terms of the characteristic dimensionless strain metric amplitude h) is h oc VIZ, where 
1Z is the merging rate of these binaries. The most "dangerous" from the point of view of 
"contamination" of the LISA frequencies are galactic binary white dwarfs with the highest 
galactic merging rate K wd ~ 1/300 yr _1 . For one- year observational run they will make up a 
stochastic GW background up to frequencies around 1 mHz (Hils and Bender 1998, Kosenko 
and Postnov 1999). 

However, if a binary star has a non-circular orbit, it emits GW at a number of harmonics 
to its orbital (Keplerian) frequency, and above estimates should be modified. This is precisely 
the case of binary neutron stars, since they generically should reside in highly eccentric orbits 
due to two supernova explosions and large mass loss in the progenitor binary system. The 
possible asymmetry of supernova explosion with a natal kick velocity of the newborn NS 
additionally affects the orbital parameters by increasing orbital eccentricities. Namely, it 
is not clear from the very beginning up to what frequency one can expect such binaries to 
form a stochastic background in one year observation (i.e. individual harmonics cannot be 
resolved within the frequency resolution bin Av ~ 3 x 10~ 8 Hz). The answer to this question 
is needed to ensure that indeed there are frequency "windows" inside the LISA frequency 
band in which a cosmological GW background can be detected. 



2 Gravitational waves in eccentric orbit 

We start with the simple case of two point masses in a circular orbit. The energy losses 
caused by the quadrapole gravitational wave emission is 

dE . . 32 G^m\m\ {m x + m 2 ) 

L = —(is K ) = —— 5 1 

at 5 & a° 

where G is Newton's gravitational constant, c is the speed of light and mi,m 2 are the masses 

of stars. All GW radiation in this case is emitted at the second harmonic to the Keplerian 

frequency. 

For a non-circular orbit with eccentricity e, the GW luminosity increases: 

" L ° x (i - e y/ 2 (2) 

The emission now is widely spread over higher-order harmonics to the Keplerian frequency 
(Peters & Mathews 1963): 

dE 

— (u = nu k ) = L x g(n, e) (3) 

n 4 2 
9( n , e) = —([J„- 2 (ne) - 2eJ n _ 1 (ne) + - J n {ne) 
32 n 

+2eJ n+1 (ne) - J„ +2 (ne)] 2 + (1 - e 2 )[J„_ 2 (ne) 

4 

-2J n (ne) + J n+2 (ne)] 2 + — J n 2 (ne)) (4) 
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Here J n is Bessel functions and n is the number of harmonics. The harmonic at which 
the maximum of energy is radiated increases with eccentricity, in accordance with the third 
Kepler's law: 

n max w (1 - e 2 )- 3/2 (5) 
This number is asymptotically precise for large eccentricities. 

3 Merging neutron star binaries 

As we noted above, an important quantity is the upper frequency vu m above which each 
individual source from a given population of binaries can be resolved during one- year obser- 
vation time (i.e. within the frequency bin Av = 3 x 1CT 8 Hz). In the case of binary white 
dwarfs, this limiting frequency is of order of ~ 1.2 x 1CT 3 Hz for galactic binary WD coa- 
lescing rate 1Z = l/300yr _1 . The coalescing rate of binary NS in most optimistic scenarios 
is about 1Z = 10~ A yr~ 1 so the limiting frequency would be vu m = 3 x 10~ 4 Hz if they were 
in circular orbits from the very beginning. However, it is not the case. 

There is observational evidence for several galactic NS binaries (observed as binary ra- 
diopulsars), all of them reciting in eccentric orbits (see the Table). The total number of NS 
binaries in the Galaxy can be roughly estimated as IZns x T ga i = 10~ 4 x 10 10 ~ 10 6 , where 
Tg a i is the galactic age. In compact binaries, the gravitational radiation back reaction is the 
only driving force for orbital evolution. Given the initial distribution of these system over 
orbital parameters, F in (a, e), it is straightforward to calculate their steady-state distribution 
(see Buitrago et al. 1994 for more detail). 



Binary PSR with NS secondaries 



PSR 


Period(days) 


eccentricitye 


t coalescing j 


B1913+16 


0.323 


0.617 


1.0 x 10 8 


B1534+12 


0.420 


0.274 


1.0 x 10 9 


B2721+llc 


0.335 


0.681 


8.0 x 10 7 


B2303+46 


12.340 


0.658 


1.6 x 10 12 



4 Delta-function examples 

Some features of interest here can be learned already from an (un realistic) delta-function- 
like initial distribution of systems: F in (a,e) = 5 (a — ao,e — eo) (hereafter we assume the 
NS+NS coalescing rate to be TZns — 10~ 4 yr _1 ). The stationary stochastic GW background 
calculated is shown in Fig. [j] and Fig. | for ao = 5-R© and eo = 0.5, 0.9, respectively. The 
contribution of first three harmonics is shown by separate curves. It is seen that increase in 
the initial eccentricity strongly affects the shape of the background up to some frequency at 
which the eccentricities of binaries in the stationary distribution become sufficiently small. 
Above this frequency only the second harmonics from almost circular binaries contributes 
to the total signal. The non-monotonic dependence of the background is due to the non- 
monotonic behavior of the energy emitted at every harmonics with eccentricity. 
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Figure 1: The GW background from a model stationary population of binary NS with initial 
delta-function like distribution at = 5R & , = 0.9 and the contribution from the first, 
second and third harmonics. For comparison, the would-be background from initially circual 
NS binaries coalescing with the same rate is shown. 




Figure 2: The same as in Fig [l| for systems with initial ao = 5R Q , eo = 0.5. 



5 Stationary stochastic GW noise from coalescing bi- 
nary NS 

Now we turn to calculation of a realistic GW background produced by binary NS. This can 
be done with the following steps. 

(1) Calculate the initial galactic distribution of binary NS stars over orbital semi-major 
axes and eccentricities F in (a, e) using e.g. binary population synthesis method (Lipunov et 
al. 1996). 

(2) Calculate a stationary distribution function F st (a,e) (Buitrago et al. 1994, Pierro V. 
and Pinto I.M 1996). 

(3) Compute the total stochastic GW background from these sources. 

We shall assume for simplicity that all the sources are at one and the same distance 
r = 7.9 kpc away, which is close to the average distance to stars in our Galaxy. This 
simplifying assumption does not change our general conclusions significantly. Basically, at 
each frequency v we sum up the GW flux from all the harmonics that fall within the chosen 
frequency bin Av = 3 x 10~ 8 Hz from lower-frequency non-circular systems in the calculated 
stationary distribution F st (y, e). 



# = TY £ I 9(n, e) F*(g, e) de) f 

u ' v \ n=l / 

n 



(6) 



£ JL Qg(n,e)F st (%,e)de) 

n=l / 

In the last equation we take into account Av <C v. The amplitude of high-order harmonics 
rapidly decreases, so we stopped the summation of harmonics for n > nu m , n lim : g(nn m , e) = 
e x max n [g(n, e)]. We assumed e = 10~ 4 , 10~ 6 , 10~ 10 . Decreasing e increases the total number 
of harmonics which contribute to the given frequency bin, but practically does not change 
the number of the most powerful harmonics within it (see Fig. |5|). 

The resulting noise curve is shown in Fig. |3| in terms of dimensionless amplitude h 

h » = -T^§ ( 7 ) 
c A r l yRv) 1 at 

As expected, the NS+NS confusion noise lies below WD+WD curve, mainly due to lower 
71ns- High-order harmonics from non-circular NS binaries mostly contribute at lower fre- 
quencies, so starting from v ~ 10~ 4 Hz the calculated noise curve practically coincides with 
that formed by circular NS binaries coalescing with the same rate TZns — 10~ 4 yr _1 . 

Extragalactic NS binaries would also form an isotropic confusion noise, but the level of 
extragalactic GW backgrounds is generally an order of magnitude lower than the galactic one 
(e.g. calculations of Kosenko and Postnov 1998) and is beyond reach by the expected LISA 
sensitivity (the bottom dashed curve in Fig. |3|). The limiting frequency vu m for extragalactic 
NS binaries finds from the condition for circular systems and can be as higher as ~ 0.3 Hz 
(Ungarelli and Vecchio 2000). 
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Figure 3: GW background from a stationary binary NS population calculated for a model 
spiral galaxy with the NS+NS coalescing rate Uns = 10~ 4 yr -1 (the middle broken curve). 
The limiting frequency above which the most powerful harmonics can be resolved in one- 
year observation is vu m ~ 0.8 mHz. If no harmonics were subtracted from the frequency bin 
3 x Hz, the background would continue toward higher frequencies as shown by the long- 
dashed line. The lower dashed line shows contribution from extragalactic NS+NS binaries. 
The upper solid line schematically represents WD+WD confusion limit from Bender and 
Hils (1997) {vum — 1.2 mHz, Eq. (1)). The proposed LISA sensitivity is from Thorne (1995) 
[Fig. 14]. 

6 The limiting frequency 

More important in the non-circular case is the increase of the limiting frequency vu m above 
which individual sources can be resolved during a one-year observation. Formally, the number 
of harmonics inside the frequency bin Av. is more than one at any frequency, and in this sense 
the GW background extends up to very high frequency ~ IKHz. But at high frequencies, 
the main contribution must come from circular binaries (2-nd harmonics), with the total 
GW power from higher-order harmonics from low-frequency eccentric systems becoming 
gradually less and less. So the number of harmonics contributing in a given frequency 
bin Av = 3 x 10~ 8 Hz were counted starting with the strongest one and continuing until 
99bin (Eq.(^)) had been included. The number of such harmonics as a function of frequency 
and the assumed neutron star kick velocity amplitude (which mostly affect the population 
synthesis results) N(v, Av) is shown in Fig. The limiting frequency vu m is determined from 
equation N(vu m , Av) = 1. To within uncertainties of our calculations (JZns, F in (a, e), etc.), 
v lim ~ 1CT 3 Hz, close to the break in the confusion limit from binary WD. For comparison, 
in the same Fig. [| we show N(u, Av) for stationary circular NS binaries coalescing with the 
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same rate TZns — 10~ 4 yr _1 (i/; jm (e = 0) oc 3 x 10~ 4 Hz). The effect of increasing the chosen 
level for the GW power from the strongest harmonics in the bin is illustrated in Fig. |5|. 
Increasing this level from 95% to 100% changes the limiting frequency vu m by almost an 
order of magnitude. 

This procedure, however, is not complete. Indeed, consider a frequency bin next to thus 
defined Pu m . We may ask the question which harmonics will be most probably absent inside 
it. The answer is those which are the least probable at this frequency. The probability to find 
a harmonics in the bin at a given frequency is determined by the number of the harmonic and 
the stationary distribution function of sources F st (i>,e). It happens that this is harmonics 
number one (mostly due to a steep decrease of the systems' stationary distribution function). 
If we remove all 1st harmonics of systems with orbital frequencies falling into the chosen bin 
from the sum Eq.(|6|), we are left again with some (smaller) GW power in the bin and may 
wish to find the limiting frequency exactly in the same way as described above (i.e. by fixing 
some level and summing up the strongest harmonics up to this level), call it fum-i- Above 
this new limiting frequency, we can repeat the entire procedure to find vum-2 (this time the 
2d harmonics is least probable to be found in the bin next to the new limiting frequency), 
etc., until the noise level of the detector is reached. The step-like line which continues the 
spectrum above ~ 1 mHz in Fig. |3] illustrates this procedure. This line of course do not 
represent the "real" GW noise curve from binary NS and just gives a feeling of how it most 
probably behaves at Vu m , above which individual harmonics from coalescing binary NS can 
be singled out. 

7 Discussion and conclusions 

A stochastic GW background is frequently described by the quantity Q,qw(v), defined as 
the ratio of the energy density GW in the bandwidth equal to frequency to the total energy 
density to close the Universe (see critical discussion of this quantity in Grishchuk et al. 
2001). In terms of this quantity, any stochastic cosmological GW background with Vt G w — 
10~ n -T- (a few) x 10~ 12 can be detected by LISA interferometer within the frequency range 
from ~ 10~ 3 to ~ 10~ 2 Hz (see Fig. (|)). 

The primordial nucleosynthesis arguments give an upper limit for possible cosmological 
GW density flow < 10~ 6 . The so-called "Standard Inflationary Model" predicts the energy 
density of relic gravitational waves at an even smaller level VLqw ~ i a f ew ) x 10~ 14 . So 
weak a background cannot be certainly detected by LISA. However, more realistic appears 
the model (Grishchuk et al. 2001) in which the cosmological stochastic gravitational wave 
background level can be well above the LISA sensitivity limit. Detecting such a background 
by the LISA interferometer remains one of the most exciting goals of this challenging project. 
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Figure 4: Number of harmonics from non-circular galactic NS binaries producing the stochas- 
tic noise in the frequency bin Au = 3 x 1CT 8 Hz as a function of the neutron star kick velocity. 
The curves are labeled with kick velocity amplitudes 100, 200 and 300 km/s assumed in the 
population synthesis calculations. The straight line is for the circular NS binaries. 
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Figure 5: Number of the most powerful harmonics from non-circular galactic NS binaries 
which contribute 95%, 99% and 100% of the total GW energy collected in the bin Au = 
3 x 10 -8 Hz. The effect of increase in the number of harmonics nu m determining the total 
energy collected in the bin is shown for different e. Kick velocity amplitude 200 km/s. 
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Figure 6: Levels of equal VLqw i n comparison with the expected stationary galactic NS+NS 
stochastic background in the LISA frequency band. 
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